Analytic Quantization of the QCD String 



Theodore J. Allen 

Physics Department, Hobart & William Smith Colleges 
Geneva, New York 14456 USA 

Charles Goebel and M. G. Olsson 

Department of Physics, University of Wisconsin, 
1150 University Avenue, Madison, Wisconsin 53706 USA 

Sinisa Veseli 
Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, Illinois 60510 USA 
(February 1, 2008) 

We perform an analytic semi-classical quantization of the straight QCD string with one end fixed and 
a massless quark on the other, in the limits of orbital and radial dominant motion. We compare our 
results to the exact numerical semi-classical quantization. We observe that the numerical semi-classical 
quantization agrees well with our exact numerical canonical quantization. 



I. INTRODUCTION 

The purpose of this paper is to explore some remark- 
able results of the QCD string/flux tube model [[l|,|J|]. 
In its relativistic and single quantized form [|^] a partic- 
ularly simple pattern emerges when one or both quarks 
are light. For the purposes of this paper we take the light 
quarks to have zero mass. In this case the energy E of the 
light degrees of freedom (LDF) and the angular and radial 
quantum numbers (J and n respectively) are accurately 
related by 



(2)7ra 



J + 2n+- 



(1.1) 



for one or (two) light quark(s). For the case of one light 
quark the meson energy is the sum of E and the heavy 
quark mass. 

The above pattern of angular and radial states results in 
degenerate "towers" of mesons of the same parity. This 
is the same pattern as the 3D harmonic oscillator. 

The QCD string model is kinematically intricate and it 
would seem u nlike ly that it would lead to such a simple 
result as Eq. (1.1). We demonstrate in this paper that 
although Eq. ( 1 . 1 ) is not exact, it is very accurate for most 



accessible quantum states. The quantized relativistic flux 
tube model is of great interest because of the probability 
that QCD reduces to string-like behavior at large source 
separations [^. 

The simplest version of a quark string model assumes 
that the string is always straight. In the limiting cases 
of circular motion or pure radial motion this assumption 



is physically reasonable. In addition, for massive quarks 
at the ends the relativistic corrections have been shown 
[^] to agree with the Wilson loop description of QCD 
confinement [^j]. Numerical quantization of the straight 
string and quark system has been done canonically 
with the Nambu-Goto string, as well as in the WKB 
approximation with an auxiliary field method [^. Both 
approaches give results similar to Eq. (1.1). 

More generally, one may allow the string to curve adia- 
batically by incorporating the string equations of motion 
from the Nambu-Goto action [|[]. This more general 
string calculation demonstrates that the string curvature 
remains small for motion in ordinary hadrons. The string 
curvature only slightly changes the energies of the bound 
states, justifying the use of the simpler straight string 
approximation. 

The approach emphasized here is to quantize the 
straight string system semi-classically. We will show 
that this quantization agrees well with our previous exact 
quantization method. We will also show that a single inte- 
gral function then predicts the whole spectroscopy when 
at least one quark is massless. We will approximate an- 
alytically this integral in sectors when J ^ n and when 
J. The later case where radial motion dominates is 
valid over most of the allowed bound states. 
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II. DYNAMICS AND QUANTIZATION 
A. Dynamics 

As we mentioned earlier, the straight QCD string is an 
excellent approximation to the dynamical (curved) string 
in normal hadrons. We will therefore restrict ourselves 
to the relatively simple straight string. We will explicitly 
consider the case of one fixed end and a quark of mass m 
at the other. The string with two light quarks introduces 
only minor modifications, which we discuss at the end 
of section ^ As is well-known [^, the two constants 
of motion are the orbital angular momentum J and the 
energy E of the light degrees of freedom, which are given 
in terms of the quark's transverse velocity v± and radial 
momentum as 



J = Wr"/±v±r + ar^f{vj_) 



(2.1) 
(2.2) 



where Wr = \/pl + rnP-, 7_l = (1 — v\)^^/-, r is the 
string length, and a is the string tension. The functions 
/ and S that appear in the expressions for the angular 
momentum and energy are 



sin^' 



(2.3) 
(2.4) 



For our present purposes we will introduce a set of 
dimensionless variables. As our units we take the circular 
orbit radius (in the limit of a massless quark) 



and corresponding string energy 

Ea = V JTTfl . 



(2.5) 



(2.6) 



Our dimensionless variables A, p, and w are defined 



by 



E 
r 



— = V1+A2 + pA , 



n) 



and 



w 



Wr 



(2.7) 
(2.8) 

(2.9) 



The leading (classical) Regge trajectory corresponds to 
A = and radial excitation occurs for positive A. 

In terms of these dimensionless variables the conserved 



quantities (2.1) and ( 2.2 ) become 



1=2 w^±_ v±_ 



Vl+A2 + pA 



W7i 



(2.10) 
(2.11) 



After some work, we can eliminate the product w^^ 
to obtain 



2 t)j_ 



Vl+A^ + pA VTTa2 



- {v±S{vA_) - f{vi_)) [Vl+A2 + pA 



= 1 



(2.12) 



and rewrite Eq. ( [2.1 1| ) to find an expression for w 

V^l +A2 



w 



S-K) ( Vl +A2 + pA 



(2.13) 



At the radial turning points pr — 0. The radial velocity 
also vanishes, except in the massless quark limit. 

The reason for the definition of the radial coordinate 
p in Eq. (2.8) will be evident in this limit. If we take 



the limit v_\_ — \ this will correspond to r = in the 
massles s qua rk case. Using /(I) — 7r/4andS'(l) = 7r/2 
in Eq. ( [2.12| ) we see that the solution is = j for any 
A. The points p = ±1 therefore correspond to possible 
turning poi nts. 

By Eq. ( 2.13 ) we see that at = 1, w = 0. This 
would seem to be the expected result in the massless 
quark case. However, as we see in Fig. |l[ the second 
factor in Eq. ( 2.13| ) has a zero at p < 1 and is the true 
turning point. In Fig. |l]we show representative numerical 
solutions with A = 0.5 and A = 2.0. First we solve 
Eq. (2.12) for v±_ for — 1 < p < +1 and substitute into 
Eq. (2.13) to compute w{p). The turning points of the 
general motion are p_ = — 1 and p+ < 1 . We will show 
later that in the massless quark case 



- 1 < P+ < 1 



(2.14) 



where the lower limit is achieved for large A and p+ = 1 
when A = 0. 

The outer turning point in the massless case is a 
"bounce" where f discontinuously changes sign. 
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B. Quantization 

The semi-classical quantization condition for a spher- 
ically symmetric system is [p|l 



dr p,. (r, J + 1/2) — tt I n + 



(2.15) 



where n = 0, 1 , 2, . . . and it is understood that the Langer 
correction replaces the classical angular momentum J 
by J + 1 /2, where J is now the angular momentum quan- 
tum number. In terms of our dimensionless parameters p 
and w, the above quantization condition becomes 



C?/9 4 / 



{2n + 1) 
(J+ i)a7r ~ (J+ i)A 



(2.16) 



For our purposes, the massless quark will be of central 
interest. In this c ase, th ere is no J dependence on the left 
hand side of Eq. (1.16) and, by defining 



4 rp+ 

/(A) = - d 

J-l 



dp w{p) , 



we obtain 



J 



(2.17) 



(2.18) 



The remarkable aspect of the massless quantization 
condition is that the whole spectrum is revealed once 
/(A) is computed. In Fig. 0, we show the numerical 
resul t of this integration. For a given A and p we use 



Eq. (2.12) to find v±. The result is use d in E q. (2.13) to 



find w{p). This is turn is used in Eq. ( 2.17 ) to compute 
/(A). The upper integration limit is determined by the 
zero in w{p), as shown in Fig. |l[ 

It is evident from Fig. ^that /(A) is essentially linear in 
A with unit slope. The LDF energy dependence is related 
to A by Eq. (O) 



£2 



{J+^)aTT 



1/2 



(2.19) 



where we have again used the Langer correction to the 
angu lar m omentum. Using the quantization condition, 
Eq. (2.18), we can map out the entire Regge structure as 
shown by the curves in Fig. ||. The trajectories are labeled 
by different radial excitations n = 0, 1,2, . . .. The n = 
trajectory is normally called the "leading trajectory" 
and the n > trajectories are known as the "daughter 
trajectories." The solid points are the exact numerical 
solutions by the exact canonical quantization method S\ . 



The nearly perfect agreement between the semi-classical 
and canonical quantization should be noted. It is only 
along the leading trajectory that small differences arise. 
If there are differences, this is where they should appear 
since the semi-classical quantization becomes exact for 
large radial excitation where the wave function has many 
nodes. 

Figure || shows that for A > 1 , / (A) is essentially equal 
to A. In general we expect an expansion in A^ ' to have 
odd powers of the form 



(2.20) 



A detai led ex amination of our exact numerical integration 
of Eq. (2.17) shows that 



-0.006 . 



(2.21) 



The small magnitude of the coefficient c reflects the ex- 
traordinary accuracy of /(A) = A for large Delta. 

The string with two light quarks can be considered 
as two single light quark strings with their fixed ends 
coinciding at the center of mass point. The resulting light 
meson equations following from Eqs. (2.1) and ( 2.2 ) are 



El ^2Wrj±+arLS{v±) , 

where = 2r, Jl — 2 J, and El = 2E. 
We recover Eq. (2.12) upon rescaling 



Eql = \/2naJL 

■A)L = 2\ ■ 
V an 



Using Eq. (2.13), we may now identify 



Wr = ^\/2naJL w{p) 



(2.22) 
(2.23) 



(2.24) 
(2.25) 

(2.26) 



and the quantization condition with /(A) ~ A yields 

Jl + 2n + 5 , (2.27) 



3. 

2'Ka 



in agreement with Eq. (1.1) 



III. RADIAL DOMINANCE ANALYTIC SOLUTION 

A. Leading Order 

As verified by examining Fig. large n corresponds to 
large A since at fixed angular momentum, both make the 
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LDF energy E increase without bound. By examining 
Eq. ( 2.12 ), we see that for large A either p = — I or v± is 
small. To verify this we note that 



v^S{v^) - /K) = 2v^/3 + 0{vi) 



(3.1) 



Hence, except at p = — 1, ^ 1 in the A ^ 1 limit. 

In the vanishing v± limit the expression (2.13) for u'(p) 
then becomes 



w[p) — > A 



1 



(3.2) 



We note that the zero in w that follows from Eq. (3.2) is 



2-' 



(3.3) 



In Fig. H we show the exact numerical result (solid 
curve) as well as the large A approximation (LO) for 
A = 2. 

In this regime we may an alytic ally evaluate the quan- 
tization integral /(A) in Eq. (2. 17) to be 



/(A)=A 



(3.4) 



Sub stituti on of this into Eq. ( |2.18| ) with the use of 
Eq. (2.19) quickly yields 



-S^ X . 3 
— ^ J + 2n+- 
ira 2 



(3.5) 



This is a good representation of the massless quark spec- 
troscopy, as we will discuss in the conclusion. Sec. ^ 
This is a reflection of the agreement of the asymptotic 
(large A) values of /(A) down to small A. 



Finally, we note that the spectrum in Eq. (3.5) is very 
similar to the three dimensional harmonic oscillator spec- 
trum as we have pointed out previously in the context of 
the Lorentz scalar confinement of massless quarks. 



B. Corrections to Leading Order 

Although in leading order in A the result 1(A) = A 
is very accurate, the distribution w{p) computed to the 
same approximation is not as satisfactory, which can be 
seen in Fig. |[ In this section we attempt to obtain a 
somewhat bet ter ap proximation to these two quantities. 

Using Eq. ( [2.12| ) and the limiting behavior of Eq. (3.1), 
we find the dependence of v± on p and A 



v± 



2 (Vl+A2 + pA 

(yrr^ - ^ ( vtta2 + pa 



(3.6) 



The value for upper turning point coordinate p+ can 
be found from the vanishing of Eq. ( 2.13 ) for small v±, 
which yields 



P+ 



yi + A^ /vr 
A V2 ^ 



(3.7) 



With this value of p+, Eq. (3.6) determines the value of 
v± at the upper turning point 



v; = I - 

TT 



1 



1 +A2 ■ 

To second order in v±, w{p) becomes 
„2 r 



(3.8) 



VT+A^- —i^vT+^ + pA 



(3.9) 



Upon substitution of the result of Eq. (3.6), becomes 



8 [VTTa2 - ^ (\/TTa2 + pA 

Wi + W2 , 



(3.10) 



with W] being the piece of Eq. ( [3.9| ) that is finite in the 
limit of vanishing v±. The next-to-leading-order result is 
shown in Fig. || We see that it comes much closer to the 
exact w{p) but deviates slightly near p = — 1 where the 
small v± expansion fails. 



It is not difficult to evaluate the integral in Eq. ( 2.17 ) 
by changing the integration variable from p to 



= vT+A^ + pA . 



(3.11) 



Splitting /(A) = /i(A) + /2(A) according to Eq. ( p.KJl ) 
we find 

2 



/i(A) 
/2(A) 



12 

A V TT 



(3.12) 



1 



27rA(l +A2) 



- - 1 - A- - aVi +A2 

TT 



4 

3^ 



37rA 



37r-4 



Asymptotically 1(A) = Ii(A) + /2(A) is Hnear 



1(A) 



1 - - 1 A- 

TT 



0(A-Hn(A)) . (3.13) 



Although the coefficient of 1 /A is small, yet higher order 
correctio ns are expected to reduce it to the near zero value 
found in (2.21 ). 
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IV. ANGULAR DOMINANCE ANALYTIC SOLUTION Numerically, we find 



For large and small radial excitation, A is small 
as can be seen from Eq. (2.18). Since for a massless 
quark the leading classical Regge trajectory corresponds 
to circular motion and in that case A = and we are 
in the ultra-relativistic regime vj « 1. In this section 



we examine the behavior of Eqs . (2.12) and ( 2 . 1 3 ) as vj 
approaches unity. 

The appropriate expansion parameter in this case is 



y = \n 



1 

-f± 



(4.1) 



which would be the radial velocity with a massless quark 
since a massless quark must hav e = + v\ — 1 . 



Our first task is to expand Eq. ( 2.12 ) for small y and A, 
which is an expansion about circular orbits. The result is 



pAy2 _ (1 _ p2)A2 = 



(4.2) 



The cubic term is critical. Without it, there would be no 
real solutions for < 1. Going to quartic order in y 
makes only minor changes. To demonstrate the accuracy 
of this appr oximation, we show in Fig. I the solution y{p) 
of Eq. ( |4.2[ ) and the exact result obtained by solving the 
original (exact) Eq. ( 2.12 ) for the value A = 0.1. 

For small A the y{p) solution becomes more symmetric 
about p = 0. The approximate solution to the cubic 



equation (4.2) for small A is 



1/3 



(4.3) 



By substitution of Eq. (4.3), we see that Eq. (4.2) is sat 



isfied up to the two leading powers of A. As A becomes 
small the asymmetric term becomes of less relative impor- 
tance and we may keep only the leading term in Eq. ( 4.3 ), 



y{p) 



A-CI 



37r 



1/3 



(1 



,1/3 



,2/3 



(4.4) 



For A <C 1 and y <C 1, the expression (2.13) for w 
reduces to 



1 



?z; ~ -y- ~ - I - ) \l~ p- 

TT 2 \TT 



2/3 



A4/3 



(4.5) 



With this approximation for w, the semi-classical quan- 
tization integral (2.17) becomes 



/(A) ~ 1.3367A'*/^ . (4.7) 

In Fig. we show the small A values of /(A). At small 
A we observe a distinct deviation from the linearity ob- 
served at larger A (see Fig. H). The dotted curve shown 
in F ig. n \ is the small A leading term approximation of 
Eq. ( |4.7f ^and the dashed curve is the large A approxima- 
tion to /(A). Although the small A analytic approxima- 
tion clearly is converging to the correct numerical result 
for A ^ 0, deviations from the exact 1(A) are evident 
around A = 0.3. For A > 0.4, the large A approximation 
to 1(A) is better 



V. CONCLUSION 

We have argued in Section | that the straight string 
model of confinement is accurate for ordinary hadron dy- 
namics and that it is therefore a serious candidate for long- 
range QCD. Semi-classical quantization of the straight 
string agrees well with a fully quantized calculation P], 
as shown in Fig. ^. The semi-classical method allows 
considerable analytic insight into string confinement in 
hadrons. In particular, when one or two of the quarks 
are massless, the entire spectroscopy is generated by a 
single integral function 1(A) given in Eq. (2.17). For 
almost all accessible states 1(A) ~ A which, by (2.18), 
immediately gives the simple pattern of Eq. (1.1), as il- 
lustrated in Fig. ||. Since even low-lying meson dynamics 
is dominated by the confining region, the observation of 
degenerate towers of mesons of the same parity becomes 
a key prediction of QCD. 



The deviations from the simple pattern (1.1) are small 
and come when A < 1 . For observed meson states this 
corresponds to large angular momentum and small ra- 
dial excitation. That is, deviations from (1.1) occur for 
mesons lying near the "leading" Regge trajectory. We 
showed in Sec. |^ that 1(A) is proportional to A"^/^ for 
A ^ 1 and the analytic approximation agrees well with 
the exact numerical 1(A), as can be seen in Fig. 0. 

To make the above statement more explicit we show 
in Fig. H a recreation of Fig. || but now adding a dashed 
line representing the analytic result (3.5) which becomes 
exact for large n and small J. For large J and small n the 
Regge trajectories seem to have a slightly smaller slope. 
In Fig. ^ we show the A values in this orbital dominant 
regime and we see that A is less than one. As seen in 
Fig. ^ 1(A) < A and this reduction accounts for the 
difference between the dashed and solid (exact) curves in 
Fig.|. 

A final observation is that 1(A) A follows from 
the radially dominant regime and not from the nearly 
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circular approximation. This limit automatically gives 
straight, evenly spaced Regge trajectories. 



ACKNOWLEDGMENTS 

This work was supported in part by the US Department 
of Energy under Contract No. DE-FG02-95ER40896. 



[1] Y. Nambu, "Quark model and the factorization of the 
Veneziano Amplitude," in Symmetries and quark mod- 
els, ed. R. Chand, (Gordon and Breach, 1970); T. Goto, 
Prog. Theor. Phys. 46, 1560 (1971) ; L. Susskind, Nuovo Ci- 
mento 69A, 457 (1970); A.M. Polyakov, Phys. Lett. 103B, 
207 (1981). 

[2] A. Chodos and C. B. Thomas, Nucl. Phys. B 72, 509 (1974); 
I. Bai-s, Nucl. Phys. B 111, 413 (1974); I. Bars and A. Hanson, 
Phys. Rev. D 13, 1744 (1976); K. Kikkawa and M. Sato, 
Phys. Rev. Lett. 38, 1309 (1977); M. Ida, Prog. Theor. Phys. 
59, 1661 (1978); K. Johnson and C. Nohl, Phys. Rev. D 19, 
291 (1979). 

[3] D. LaCourse and M. G. Olsson, Phys. Rev. D 39, 2751 

(1989) ; C. Olson, M. G. Olsson, and K. WiUiams, Phys. Rev. 
D 45, 4307 (1992); M. G. Olsson and K. Williams, Phys. Rev. 
D 48, 417 (1993); M. G. Olsson and S. Veseh, Phys. Rev. D 
51, 3578 (1995); M. G. Olsson, S. Veseli, and K. Williams, 
Phys. Rev. D 53, 4006 (1996). 

[4] T.J. Allen, M.G. Olsson, and S. Veseli, Phys. Lett. B 434 

(1998) 110. 

[5] C. Olson, M.G. Olsson, and K. Williams, Phys. Rev. D 45, 
4307 (1992). 

[6] N. Brambilla and B.M. Prosperi, Phys. Lett. B 236, 69 

(1990) ; G.S. Bah and P Boyle, Phys. Rev. D 59, 114504 

(1999) . 

[7] A.Yu. Dubin, A.B. Kaidalov, and Yu.A. Simonov, 
Phys. Lett. B 323, 41 (1994); Yad. Fiz. 56, 213 (1993); 
V.L. Morgunov, A.V. Nefediev and Yu.A. Simonov, Phys. 
Lett. B 459, 653 (1999). 

[8] T.J. Allen, M.G. Olsson, and S. Veseli, Phys. Rev. D 59 
094011 (1999); T.J. Allen, M.G. Olsson, and S. Veseli, 
Phys. Rev. D 60 074026 (1999) 

[9] R.E. Langer, Phys. Rev. 51, 669 (1937); M.S. Child, Molec- 
ular Collision Theory (Academic, New York, 1974). 
[10] T.J. Allen, M.G. Olsson, and S. Veseli, Phys. Rev. D 62 
094021 (2000). 



Q. 



2.0 
1.5 
1.0 
0.5 
0.0 
■0.5 
■1.0 



y \x A=2.0 






A = 0.5 


A = 0.5 ^ 


\ A = 2.0 




w(p) 




vjp) 





-1.0 



■0.5 



0.0 
P 



0.5 



1.0 



FI G. 1 . Exact numerical solutions of Eq. ( 2.12 ) and 
Eq. ( [2.I3I ) for wIp) and v±{p) for the values A = 0.5 and 
A = 2.0. 
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FIG. 4. The dependence of A upon the radial quantum num- 
ber n and the angular momentum quantum number J. Increas- 
ing n leads to increasing A. 



FIG. 5. A comparison of the exact numerical solution of 



Eq. (2.13) for w{p) (lowest curve) with the leading order solu- 
tion Eq. (3.2) (highest curve) and the solution to next to leading 
order, Eq. (3.1G) (intermediate curve). 
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FIG. 6. A comparison of the exac t numerical solution for 
y{p) ~ (1 ~ '''i)'^" from Eq. ( |2.13| ) and the solution of the 
cubic approximation, Eq. Q,forA = 0.1. 



FIG. 7. A comparison of the exact numerical solution for 
/(A) (solid line) and the small A approximation (dotted line) 
and the large A approximation (dashed line). 
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FIG. 8. The Regge trajectories from the large A approxi- 
mation (dashed lines) of Eq. (1.1) in comparison to the exact 
numerical semi-classical analysis (lines) and the results of nu- 
merical canonical quantization (dots). 



FIG. 9. The dependence of A upon the angular momentum 
quantum number J and the radial quantum number n. Increas- 
ing J leads to decreasing A. 
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